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Abstract 

Due to the observation of Bose-Einstein condensation in dilute atomic ^^Rb 
and ^^Na vapors last year, there is currently a great interest in the properties 
of a degenerate Bose gas. One aspect that is not yet understood is the effect 
of the interatomic interactions on the critical temperature of the phase transi- 
tion. We present here a renormalization group study of the weakly-interacting 
Bose gas and predict the critical temperature at which Bose-Einstein con- 
densation occurs as a function of the (positive) scattering length a of the 
interatomic potential. 
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After a long history of trying to master stabilization and cooling, the correct road to- 
wards Bose-Einstein condensation in a dilute gas has now finally been found. Indeed, last 
year a Bose condensate was observed in the magnetically trapped and evaporatively cooled 
alkali gases ^^Rb and ^^Na Initiated by this breaktrough, extensive experimental and 

theoretical studies will surely follow in the near future to fathom the physics of Bose-Einstein 
condensation in inhomogeneous gases. Superfluidity P], the Josephson effect and vortex dy- 
namics are just a few examples of phenomena that are very interesting, but in first instance 
also very difficult to experimentally measure or create. Aspects that appear to be more 
easily accessible and will be subject of scrutiny first, are e.g. collective excitations IQ and 
the critical temperature of the phase transition. With respect to this latter point, the recent 
experiments claim to observe a critical temperature that is shifted upward compared to the 
ideal gas value, as the measured degeneracy parameter ncA^^ at the transition is smaller 
than C(3/2) — 2.612. The exact magnitude of the shift however, is still uncertain. Here n is 
the density and A^/j = (27[h'^ / mkBTy^"^ the thermal de Broglie wavelength of the atoms in 
the gas. 

A shift in the critical temperature can be expected because the gas contains a finite 
number of particles. However, this effect is negligible for the large number of atoms used 
in the experiments f^, and would shift the critical temperature downwards. In addition, 
a shift can arise because the atoms are interacting via an (effectively) repulsive potential. 
Concerning this latter aspect, it is important to note that existing approaches to the weakly- 
interacting Bose gas using the Bogoliubov (or Popov) theory actually do not predict a shift 
in the critical temperature . Since these theories are of mean-field type, one might naively 
expect that critical fluctuations will depress the critical temperature, which would again 
contradict the experimental data. To settle this issue, and also give a quantative prediction 
for the shift of the critical temperature, we report here on a renormalization group (RG) 
calculation for the homogeneous Bose gas. 

In most applications, the RG method is used to obtain information about universal prop- 
erties of the system considered, such as critical exponents. Nonuniversal properties depend 
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on the (usually unknown) ultraviolet cutoff A of the theory. However, in our case we can 
eliminate this cutoff dependence because we have sufficient information on the microscopic 
details of the dilute Bose system. Therefore, we are also able to determine nonuniversal 
properties of the system as we will see shortly. 

In principle, we treat here only the homogeneous Bose gas with effectively repulsive 
interactions, which implies that the interatomic potential has a positive scattering length 
a. However, as the number of particles ^ 1, the critical temperature Tc is much larger 
than the energy splitting huj between the levels in the harmonic oscillator traps used ex- 
perimentally. As the width of the region in which finite-size effects are important is of 
0{Tc{huj / ksTc)) , one can practically for all temperatures use a local density approximation 
to describe the gas in the center of the trap. As a result, the critical conditions we find here 
for the homogeneous gas are essentially also valid for the central density of a trapped Bose 
gas, and in particular pertain also to the ^''Rb and ^^Na experiments. 

In the functional formulation of the grand canonical partition function, the action for 
the dilute Bose system in momentum and frequency space is written as 

S[a*, a] = Y^{-ihuJn + Ck - At)ak_„ak,n 

k,n 

+ 2 f^gy X! ^q'^k+q,n+m'^k'-q,n'-m'^k',n'«k,n , (l) 

k,k',q 

n,n' ,m 

where ojn = 27m/^/3 = 2TinkBT /h are the bosonic Matsubara frequencies, ek = ^^k^/2m, n 
is the chemical potential, V is the volume of the system, V"q = / (ix y(x)e~*'^'^ is the Fourier 
transform of the interaction potential, and aj^^ and the Fourier components of the 

fields '?/'*(x, r) and ■?/'(x, r), creating respectively annihilating a particle at position x and 
(imaginary) time r. The specific form of the potential V^(x) or its corresponding scattering 
length is in first instance immaterial to the derivation of the RG equations. 

The renormalization of the various terms in the effective action is obtained by succes- 
sively integrating out high momentum shells dA of infinitesimal width 0. In the operator 
language this is equivalent to performing the trace in the grand canonical partition function 



over the most rapidly oscillating one-particle states, and then constructing a new effective 
Hamiltonian for the less rapidly oscillating states. Technically, this first step of Wilson's RG 
transformation boils down to calculating the Feynman diagrams renormalizing the vertices of 
interest, where the integration over the internal momenta is restricted to the afore-mentioned 
momentum shell. The type of diagrams contributing to the renormalization when the thick- 
ness of the momentum shell dA is infinitesimal, are the so-called one-loop diagrams only. 
This follows from the fact that each extra loop introduces an extra factor of dA The sec- 
ond step in Wilson's RG transformation consists of a rescaling of the momenta, frequencies 
and fields such that the original cutoff is restored and the coefficient of {—ihun + €k)a^ „ak,n 
is kept equal to one. If one neglects the renormalizations from the first step, this procedure 
yields the so-called trivial scaling of the vertices and reveals the relevance of the vertices 
considered. 

In principle, the action in Eq. (1) does not possess a real sharp ultraviolet cutoff A. 
However, the typical behavior of the Fourier transform of the two-body interaction potential 
is such that there is an effective ultraviolet cutoff around the momentum scale set by the 
scattering length a of this potential. Below this value, the Fourier transform is practically 
momentum independent and equal to Vq. As in the Bose systems considered here we always 
have that h/a ^ h/Ath, the momentum range in which the particles in the gas reside 
falls well below this cutoff. Thus, we represent the interaction potential by the momentum 
independent value Vq for q less than a cutoff A of 0(1 /a) and zero for g > A. 

Modelling the potential as described, implies that the nonuniversal properties we find 
from an RG calculation will in first instance be sensitive to the specific value of the cutoff 
A taken in the calculations. However, at this point our knowledge about the microscopic 
details of the Bose gas comes in to resolve this problem because we know that the two-body 
interaction potential Vq has to renormalize to the two-body T— matrix T^'^((k — k')/2 + 
q, (k — k')/2; fi^(k — k')^/m) when we include all possible two-body scattering processes in 
the vacuum. The two-body T— matrix has roughly the same momentum dependence as V^, 
and is in particular constant and equal to Airah^ /m in the range of thermal momenta and 
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energies. Thus, given an ultraviolet cutoff A, we can fix the RG equations by demanding that 
for the two-body problem Vq indeed correctly renormalizes to 'inak'^/m. Since this value is, 
due to the inequality a/Ath <C 1, already attained before we enter into the thermal regime as 
we integrate out more and more momentum shells, this indeed leads to a correct description 
of the properties of the Bose gas which is independent of the ultraviolet cutoff A. Having 
eliminated the cutoff dependence in this manner, we are then in a position to determine 
also the nonuniversal properties of the dilute Bose gas such as most notably the critical 
temperature. Furthermore, we can perform the calculation for any (positive) value of the 
scattering length, thus being able to describe any atomic species with effectively repulsive 
s— wave scattering. The results we find are therefore relevant to the experiments using ^^Rb, 
^^Na, but for instance also for future experiments with and other atoms having a positive 
scattering length. 

When deriving the equations governing the renormalization of the various terms in the 
action, one has to proceed differently, depending on the sign of the chemical potential. For 
negative chemical potential the derivation of the RG equations is most straightforward as 
one can start directly from Eq. (1) to obtain the Feynman diagrams. For positive chemical 
potential the derivation is much more involved since in this case the space and time inde- 
pendent part of the effective action, i.e. — | ao,o P +Vo | ao,o I'' /2hpV, has a Mexican hat 
shape and one first has to expand the action around the correct extremum by performing 
the shift ao,o ao,o + y/nohfiV introducing the condensate density Uq. Only after that can 
one find the Feynman diagrams contributing to the renormalization. 

We will first focus on the situation of negative chemical potential, and restrict ourselves to 
operators which are relevant or marginal at the critical temperature. The chemical potential 
jj, and the two-body interaction term Vq are both relevant there. The vertex Uq from the 
three-body interaction term, schematically written ets Uq \ t) \^ /6 and in principle also 
present in the Hamiltonian of the system, is marginal, and would thus have to be included 
in the set of RG equations describing the change or 'flow' of the vertices during the process 
of integrating out momentum shells. However, since we are studying the dilute Bose gas it it 
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expected that three-body interactions will in general have only an extremely small influence 
on the properties to be calculated. Only near the critical temperature can there be some 
effect as at that point all fluctuations become important. This we have checked explicitly by 
numerically solving the coupled set of RG equations including and excluding the three-body 
interaction. Indeed, the above argument is corroborated by our findings, which show that 
the influence of Uo on e.g. the density, even in the critical region, is typically smaller than 
1% in the regime where na^ -C 1. In the following we will, therefore, neglect the three- 
body interaction and focus on the renormalization of n and Vq only. For negative chemical 
potential, the RG equations are then given by 

^ = 2/i-^VoiV(eA-/i) (2a) 

dyo^_y_ A^' 2 l + 2iV(eA-/i) 
dl ° 27r2 ° 2(eA-/i) 

Vi/3N{eA - fi) [N{eA - ^) + 1] . (2b) 
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At a particular value of /, the region in momentum space that has been integrated out is 
the shell running from | k |= Ae~' to | k |= A. Here, A^(eA — jj) = l/(e^*^'^'^^^'' — 1) is the 
Bose- Einstein distribution function containing eA = fi^ Is?' 127X1 in its argument since after 
each step we have rescaled the momenta to retain the cutoff at A. Furthermore, as a result 
of this rescaling we have that the temperature scales as T{1) = Te^K The nontrivial term in 
Eq. (2a) originates from the sum of the Hartree and Fock diagrams, whereas the nontrivial 
terms in Eq. (2b) stem from the ladderdiagram and the bubblediagrams respectively. 

Next we consider the case of positive chemical potential. Introducing the condensate 
density through the shift mentioned above, we generate, next to a linear term proportional 
to (ao,o + c^oo)' normal and anomalous selfenergies Z]k,n(^^ii'^kn'^k,n + ^Si2('^kn'^-k -n + 

«k,na-k,-n)/2) and the term Z]k,nI]q,mr3(aq,m'^k-q,n-m'^k,n + Ctk+q,n+m'^q,»nak,n)/ VTl^V luto 

the action 0. In first instance we have, due to the neglect of the three-body interactions, 
that hT,ii = 2noVo, = uqVo, and Fs = ^/noVo■ Moreover, the magnitude of the 

condensate is found by eliminating the linear term in the action, which gives uq = fi/Vo- 
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Performing a Bogoliubov transformation to diagonalize the quadratic part of the action 
facihtates the calculation of the one-loop diagrams contributing to the renormalization of the 
vertices of interest. One of these corresponds to the linear term in the action. After every 
step of the RG transformation we have to eliminate this term by performing a shift in ao,o in 
order to remain in the minimum of the action. Clearly this gives the equation describing the 
change of the condensate density. Proceeding with the selfenergies, we find that they obey 
the Hugenholtz-Pines relation = — /iSi2(/) implicate from the U{1) symmetry 

of the problem P]. As a result, we only have to consider the diagrams contributing to /iSi2 
to obtain the renormalization of the quadratic part of the action. Finally, due to the neglect 



of the three-body interactions, we have again from U{1) symmetry that T3{1) = y no{l)Vo{l) 
and Vo(/) = hJ]i2{l)/no{l). These relations allow the action to be recast into the explicitly 
f/(l) symmetric form of Eq. (1) at any point during renormalization. Consequently, from 
the RG equations for fiT^u and no, we can deduce the equations describing the flow of the 
chemical potential fi and the two-body interaction Vq. This leads to 

dl ^ 27r2 



2ei + 6/.ei + /.3^2Ar(ac.,) + l)-l 



/i( 2eA + /i) 
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+ .2 2 f3N{huJA)[N{huJA) + 1] 



(3a) 



dVo_ 2 



\2 



where the Bose-Einstein distribution function now contains the Bogoliubov dispersion hu 



+ ' X 2 f3N{hujA)[N{huA) + 1] 



(3b) 



e\ + 2/ieA. Eqs. (2) and (3) constitute the flow equations in the case of negative and 
positive chemical potential respectively, and completely descibe the system in the dilute gas 
approximation. Note that both sets coincide when ^ is taken equal to zero, implying that 
the flow is continuously differentiable, also at yU = 0. 

Finally, we have to fix the RG equations as explained above. From Eq. (2b) we recognize 
that in a vacuum, i.e. A^(eA — yu) = 0, the renormalization of the interaction between two 
particles is governed by 



dVo 



dl 




2(eA - ii) 



1 



(4) 



As expected, this is just the differential form of the Lippmann-Schwinger equation for the 
two-body T— matrix at energy 2/x, T^^ {0,0; 2/1). As the two-body T— matrix is energy 



requirement is now that, given an ultraviolet cutoff A, Vq flows to the value Anah'^/m. This 
can be ascertained by choosing the right initial condition for Vq. More precisely we find from 
Eq. (4) that Vo{l = 0) = Anah'^ / {m{l — 2aA/7r)) leads to the correct result. Clearly, we can 
describe different atomic species by only changing the value of the scattering length a used in 
this equation. Note again that this boundary value also assures the correct renormalization 
of the full RG equations (2)- (3), because the other terms in the right-hand side of these 
equations are insensitive to the cutoff. Put differently, these terms only start to contribute 
to the flow in the range of thermal momenta, but at that point Vq has already correctly 
renormalized to A'KaU'^ /m. 

Having found the correct boundary conditions, jiil = 0) being the bare chemical potential 
in Eq. (1), we can now numerically integrate the RG equations. In Fig. 1 we depict the 
trajectories resulting from these calculations. The temperature is fixed and the bare chemical 
potential is raised going from a lower curve to a higher curve. We can determine the 
corresponding density by noting that n = no + Z]k,n('^k,n'^k,n), which we can also cast into 
a differential equation yielding the building up of the density as we integrate out more and 
more momentum shells. The trajectory flowing into the fixed point (/i*,V^*) corresponds 
after removal of the trivial scaling to a chemical potential renormalizing exactly to zero, 
and therefore to a condensate density renormalizing exactly to zero. For values higher than 
this (positive) bare critical chemical potential we find Uq > 0, and we are in the condensed 
phase. 

Thus, the critical density of the dilute Bose system is determined by the set Eq. (3). 
Moreover, we can find the usual critical exponent v pertaining to the divergence of the 
correlation length as we approach the critical temperature, i.e. ^ behaves as .^o | {T — 



independent for low energies, we can neglect fi and use T^^(0,0;2//) = Anah^/m. The 
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Tc)/Tc \~^, by linearizing the flow equations around the flxed point and calculating the 
largest eigenvalue. This gives u = 0.685, which is to be compared with the value u = 0.67 
found from an e— expansion of the 0(2) model [|T^ and also measured in ^He experiments. 
The agreement is surprisingly good, indicating that we are indeed accurately describing the 
Bose gas with the derived RG equations, also in the critical region. 

We now turn to the ultimate goal of this study, namely the effect of interactions on 
the critical temperature of Bose-Einstein condensation. At flxed temperatures we vary the 
(bare) chemical potential to flnd the trajectories flowing into the flxed point. As mentioned 
above, this yields the critical densities for these speciflc temperatures and gives us the ric — T 
relation. We repeat this for different values of the scattering length to obtain the dependence 
of the critical temperature on the strength of the interaction. All curves fall on one single 
curve if we plot the degeneracy parameter n^A^^^ at which Bose-Einstein condensation occurs 
versus a/Ath- This curve is depicted in Fig. 2, where we normalized ncA^^ to the ideal gas 
value C(3/2). As seen from Fig. 2 we always have that n^A^^ < ^(3/2), which leads to the 
conclusion that the critical temperature of the weakly-interacting Bose gas is raised with 
respect to the ideal gas value, in qualitative agreement with the recent experiments. The 
main reason for the lower value of ^cA^/j is that at the critical point, the effective chemical 
potential renormalizes from a positive initial value exactly to zero. Due to the Bogoliubov 
dispersion, this behavior depresses the occupation of the non-zero momentum states relative 
to the ideal gas case. 

In summary, we have derived renormalization group equations for the dilute Bose gas. 
Using our knowledge of the two-body scattering problem, we have effectively eliminated the 
ultraviolet cutoff from our theory, and obtained also information on nonuniversal properties. 
As such, we have determined the effect of repulsive interactions on the critical temperature 
of Bose-Einstein condensation, and found that it is raised with respect to the ideal gas value. 



Preliminary Quantum Monte Carlo calculations seem to conflrm this result |Tl|. From our 
calculations we predict that for the '^''Rb and ^^Na experiments the critical temperature can 
be raised with as much as 10 %. This appears to be a very promising result because one 



might expect that an effect of this magnitude can very well be measured in future, more 
accurate, experiments. 

We acknowledge helpful discussions with Eric Cornell, Wolfgang Ketterle and Steve 
Girvin. 
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FIGURES 

FIG. 1. Flowdiagram resulting from the RG equations. The fixed point is indicated with an 

asterix and its eigenvectors by the dotted Unes. 



FIG. 2. The normahzed critical degeneracy parameter versus a/Aj/j. 
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